function dbasis = chebyshev_dbasis(N, t)
% 计算切比雪夫多项式基函数的导数  
% 从下标0开始基函数  
% 输入参数：  
% N: 基函数的个数  
% t: 时间变量 [-1 1]  
% 输出参数：  
% dbasis: 切比雪夫多项式基函数的导数  
error("bug");
    dbasis = zeros(length(t), N);
    
    for k = 1:N
        n = k - 1;  % 切比雪夫多项式的阶数，从 0 开始
        if n == 0
            dbasis(:, k) = 0;  % T_0(x) = 1，其导数为 0
        elseif n == 1
            dbasis(:, k) = ones(size(t));  % T_1(x) = x，其导数为 1
        else
            dbasis(:, k) = n * sin(n * acos(t)) ./ sqrt(1 - t.^2);
        end
    end
end
